Multilevel Deconstruction of the In Vivo Behavior of Looped DNA-Protein Complexes

Protein-DNA complexes with loops play a fundamental role in a wide variety of cellular processes, ranging from the regulation of DNA transcription to telomere maintenance. As ubiquitous as they are, their precise in vivo properties and their integration into the cellular function still remain largely unexplored. Here, we present a multilevel approach that efficiently connects in both directions molecular properties with cell physiology and use it to characterize the molecular properties of the looped DNA-lac repressor complex while functioning in vivo. The properties we uncover include the presence of two representative conformations of the complex, the stabilization of one conformation by DNA architectural proteins, and precise values of the underlying twisting elastic constants and bending free energies. Incorporation of all this molecular information into gene-regulation models reveals an unprecedented versatility of looped DNA-protein complexes at shaping the properties of gene expression.Comment: Open Access article available at http://www.plosone.org/article/fetchArticle.action?articleURI=info%3Adoi%2F10.1371%2Fjournal.pone.000035

The transition between stochastic and deterministic behavior in an excitable gene circuit

We explore the connection between a stochastic simulation model and an ordinary differential equations (ODEs) model of the dynamics of an excitable gene circuit that exhibits noise-induced oscillations. Near a bifurcation point in the ODE model, the stochastic simulation model yields behavior dramatically different from that predicted by the ODE model. We analyze how that behavior depends on the gene copy number and find very slow convergence to the large number limit near the bifurcation point. The implications for understanding the dynamics of gene circuits and other birth-death dynamical systems with small numbers of constituents are discussed.Comment: PLoS ONE: Research Article, published 11 Apr 201

Spatial and topological organization of DNA chains induced by gene co-localization

Transcriptional activity has been shown to relate to the organization of chromosomes in the eukaryotic nucleus and in the bacterial nucleoid. In particular, highly transcribed genes, RNA polymerases and transcription factors gather into discrete spatial foci called transcription factories. However, the mechanisms underlying the formation of these foci and the resulting topological order of the chromosome remain to be elucidated. Here we consider a thermodynamic framework based on a worm-like chain model of chromosomes where sparse designated sites along the DNA are able to interact whenever they are spatially close-by. This is motivated by recurrent evidence that there exists physical interactions between genes that operate together. Three important results come out of this simple framework. First, the resulting formation of transcription foci can be viewed as a micro-phase separation of the interacting sites from the rest of the DNA. In this respect, a thermodynamic analysis suggests transcription factors to be appropriate candidates for mediating the physical interactions between genes. Next, numerical simulations of the polymer reveal a rich variety of phases that are associated with different topological orderings, each providing a way to increase the local concentrations of the interacting sites. Finally, the numerical results show that both one-dimensional clustering and periodic location of the binding sites along the DNA, which have been observed in several organisms, make the spatial co-localization of multiple families of genes particularly efficient.Comment: Figures and Supplementary Material freely available on http://dx.doi.org/10.1371/journal.pcbi.100067

Noisy-threshold control of cell death

<p>Abstract</p> <p>Background</p> <p>Cellular responses to death-promoting stimuli typically proceed through a differentiated multistage process, involving a lag phase, extensive death, and potential adaptation. Deregulation of this chain of events is at the root of many diseases. Improper adaptation is particularly important because it allows cell sub-populations to survive even in the continuous presence of death conditions, which results, among others, in the eventual failure of many targeted anticancer therapies.</p> <p>Results</p> <p>Here, I show that these typical responses arise naturally from the interplay of intracellular variability with a threshold-based control mechanism that detects cellular changes in addition to just the cellular state itself. Implementation of this mechanism in a quantitative model for T-cell apoptosis, a prototypical example of programmed cell death, captures with exceptional accuracy experimental observations for different expression levels of the oncogene Bcl-x_Land directly links adaptation with noise in an ATP threshold below which cells die.</p> <p>Conclusions</p> <p>These results indicate that oncogenes like Bcl-x_L, besides regulating absolute death values, can have a novel role as active controllers of cell-cell variability and the extent of adaptation.</p

Bistability, Probability Transition Rate and First-Passage Time in an Autoactivating Positive-Feedback Loop

A hallmark of positive-feedback regulation is bistability, which gives rise to distinct cellular states with high and low expression levels, and that stochasticity in gene expression can cause random transitions between two states, yielding bimodal population distribution (Kaern et al., 2005, Nat Rev Genet 6: 451-464). In this paper, the probability transition rate and first-passage time in an autoactivating positive-feedback loop with bistability are investigated, where the gene expression is assumed to be disturbed by both additive and multiplicative external noises, the bimodality in the stochastic gene expression is due to the bistability, and the bistability determines that the potential of the Fokker-Planck equation has two potential wells. Our main goal is to illustrate how the probability transition rate and first-passage time are affected by the maximum transcriptional rate, the intensities of additive and multiplicative noises, and the correlation of additive and multiplicative noises. Our main results show that (i) the increase of the maximum transcription rate will be useful for maintaining a high gene expression level; (ii) the probability transition rate from one potential well to the other one will increase with the increase of the intensity of additive noise; (iii) the increase of multiplicative noise strength will increase the amount of probability in the left potential well; and (iv) positive (or negative) cross-correlation between additive and multiplicative noises will increase the amount of probability in the left (or right) potential well

Exploring hypotheses of the actions of TGF-beta 1 in epidermal wound healing using a 3D computational multiscale model of the human epidermis

In vivo and in vitro studies give a paradoxical picture of the actions of the key regulatory factor TGF-beta 1 in epidermal wound healing with it stimulating migration of keratinocytes but also inhibiting their proliferation. To try to reconcile these into an easily visualized 3D model of wound healing amenable for experimentation by cell biologists, a multiscale model of the formation of a 3D skin epithelium was established with TGF-beta 1 literature-derived rule sets and equations embedded within it. At the cellular level, an agent-based bottom-up model that focuses on individual interacting units ( keratinocytes) was used. This was based on literature-derived rules governing keratinocyte behavior and keratinocyte/ECM interactions. The selection of these rule sets is described in detail in this paper. The agent-based model was then linked with a subcellular model of TGF-beta 1 production and its action on keratinocytes simulated with a complex pathway simulator. This multiscale model can be run at a cellular level only or at a combined cellular/subcellular level. It was then initially challenged ( by wounding) to investigate the behavior of keratinocytes in wound healing at the cellular level. To investigate the possible actions of TGF-beta 1, several hypotheses were then explored by deliberately manipulating some of these rule sets at subcellular levels. This exercise readily eliminated some hypotheses and identified a sequence of spatial-temporal actions of TGF-beta 1 for normal successful wound healing in an easy-to-follow 3D model. We suggest this multiscale model offers a valuable, easy-to-visualize aid to our understanding of the actions of this key regulator in wound healing, and provides a model that can now be used to explore pathologies of wound healing

Analysis of In-Vivo LacR-Mediated Gene Repression Based on the Mechanics of DNA Looping

Interactions of E. coli lac repressor (LacR) with a pair of operator sites on the same DNA molecule can lead to the formation of looped nucleoprotein complexes both in vitro and in vivo. As a major paradigm for loop-mediated gene regulation, parameters such as operator affinity and spacing, repressor concentration, and DNA bending induced by specific or non-specific DNA-binding proteins (e.g., HU), have been examined extensively. However, a complete and rigorous model that integrates all of these aspects in a systematic and quantitative treatment of experimental data has not been available. Applying our recent statistical-mechanical theory for DNA looping, we calculated repression as a function of operator spacing (58–156 bp) from first principles and obtained excellent agreement with independent sets of in-vivo data. The results suggest that a linear extended, as opposed to a closed v-shaped, LacR conformation is the dominant form of the tetramer in vivo. Moreover, loop-mediated repression in wild-type E. coli strains is facilitated by decreased DNA rigidity and high levels of flexibility in the LacR tetramer. In contrast, repression data for strains lacking HU gave a near-normal value of the DNA persistence length. These findings underscore the importance of both protein conformation and elasticity in the formation of small DNA loops widely observed in vivo, and demonstrate the utility of quantitatively analyzing gene regulation based on the mechanics of nucleoprotein complexes

Stochastic De-repression of Rhodopsins in Single Photoreceptors of the Fly Retina

The photoreceptors of the Drosophila compound eye are a classical model for studying cell fate specification. Photoreceptors (PRs) are organized in bundles of eight cells with two major types – inner PRs involved in color vision and outer PRs involved in motion detection. In wild type flies, most PRs express a single type of Rhodopsin (Rh): inner PRs express either Rh3, Rh4, Rh5 or Rh6 and outer PRs express Rh1. In outer PRs, the K50 homeodomain protein Dve is a key repressor that acts to ensure exclusive Rh expression. Loss of Dve results in de-repression of Rhodopsins in outer PRs, and leads to a wide distribution of expression levels. To quantify these effects, we introduce an automated image analysis method to measure Rhodopsin levels at the single cell level in 3D confocal stacks. Our sensitive methodology reveals cell-specific differences in Rhodopsin distributions among the outer PRs, observed over a developmental time course. We show that Rhodopsin distributions are consistent with a two-state model of gene expression, in which cells can be in either high or basal states of Rhodopsin production. Our model identifies a significant role of post-transcriptional regulation in establishing the two distinct states. The timescale for interconversion between basal and high states is shown to be on the order of days. Our results indicate that even in the absence of Dve, the Rhodopsin regulatory network can maintain highly stable states. We propose that the role of Dve in outer PRs is to buffer against rare fluctuations in this network

Persistent Cell-Autonomous Circadian Oscillations in Fibroblasts Revealed by Six-Week Single-Cell Imaging of PER2::LUC Bioluminescence

Biological oscillators naturally exhibit stochastic fluctuations in period and amplitude due to the random nature of molecular reactions. Accurately measuring the precision of noisy oscillators and the heterogeneity in period and strength of rhythmicity across a population of cells requires single-cell recordings of sufficient length to fully represent the variability of oscillations. We found persistent, independent circadian oscillations of clock gene expression in 6-week-long bioluminescence recordings of 80 primary fibroblast cells dissociated from PER2::LUC mice and kept in vitro for 6 months. Due to the stochastic nature of rhythmicity, the proportion of cells appearing rhythmic increases with the length of interval examined, with 100% of cells found to be rhythmic when using 3-week windows. Mean period and amplitude are remarkably stable throughout the 6-week recordings, with precision improving over time. For individual cells, precision of period and amplitude are correlated with cell size and rhythm amplitude, but not with period, and period exhibits much less cycle-to-cycle variability (CV 7.3%) than does amplitude (CV 37%). The time series are long enough to distinguish stochastic fluctuations within each cell from differences among cells, and we conclude that the cells do exhibit significant heterogeneity in period and strength of rhythmicity, which we measure using a novel statistical metric. Furthermore, stochastic modeling suggests that these single-cell clocks operate near a Hopf bifurcation, such that intrinsic noise enhances the oscillations by minimizing period variability and sustaining amplitude

‘Glocal’ Robustness Analysis and Model Discrimination for Circadian Oscillators

To characterize the behavior and robustness of cellular circuits with many unknown parameters is a major challenge for systems biology. Its difficulty rises exponentially with the number of circuit components. We here propose a novel analysis method to meet this challenge. Our method identifies the region of a high-dimensional parameter space where a circuit displays an experimentally observed behavior. It does so via a Monte Carlo approach guided by principal component analysis, in order to allow efficient sampling of this space. This ‘global’ analysis is then supplemented by a ‘local’ analysis, in which circuit robustness is determined for each of the thousands of parameter sets sampled in the global analysis. We apply this method to two prominent, recent models of the cyanobacterial circadian oscillator, an autocatalytic model, and a model centered on consecutive phosphorylation at two sites of the KaiC protein, a key circadian regulator. For these models, we find that the two-sites architecture is much more robust than the autocatalytic one, both globally and locally, based on five different quantifiers of robustness, including robustness to parameter perturbations and to molecular noise. Our ‘glocal’ combination of global and local analyses can also identify key causes of high or low robustness. In doing so, our approach helps to unravel the architectural origin of robust circuit behavior. Complementarily, identifying fragile aspects of system behavior can aid in designing perturbation experiments that may discriminate between competing mechanisms and different parameter sets